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The role of asymmetry on the evolution of prebiotic homochirality is investigated in the context 
of autocatalytic polymerization reaction networks. A model featuring enantiometric cross-inhibition 
and chiral bias is used to study the diffusion equations controlling the spatiotemporal development 
of left and right-handed domains. Bounds on the chiral bias are obtained based on present-day 
constraints on the emergence of life on early Earth. The viability of biasing mechanisms such as 
weak neutral currents and circularly polarized UV light is discussed. The results can be applied to 
any hypothetical planetary platform. 
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O ■ I. INTRODUCTION 

<N : 

[ The emergence of biomolecular homochirality in prebiotic Earth is a crucial step in the early history of life [J [B] ■ It 
is well-known that chiral selectivity plays a key role in the biochemistry of living systems: amino acids in proteins are 
left-handed while sugars are right-handed. However, laboratory syntheses produce racemic results. This is somewhat 

■ surprising, given that statistical fluctuations of reactants will invariably bias one enantiometer over the other even 
T— I ' though every synthesis is ab initio asymmetric 0, the enantiometric excess is nevertheless erased as the reactions 

unfold. An important exception is the reaction by Soai and coworkers, where a small initial enantiometric excess 

■ is effectively amplified in the autocatalytic alkylation of pyrimidyl aldehydes with dialkylzincs [31} . As stressed by 
, Blackmond Soai's reaction succeeds because it features the needed autocatalytic behavior proposed originally by 
' Frank with enantiometric cross-inhibition catalysed by dimers. 

ly-j \ It is unlikely that the specific chemistry of the Soai reaction occurred in early-Earth. However, it displays the relevant 
\^ • signatures of a realistic homochirality-inducing reaction network: autocatalysis, enantiometric cross-inhibition, and 

■ enzymatic enhancement performed by dimers or by larger chirally-pure chains. In the present work, we will investigate 
. the spatiotemporal dynamics of a reaction network recently proposed by Sandars which shares these features [30] . To 
' it, we will add an explicit chiral bias, in order to investigate the efficacy of intrinsic and extrinsic biasing mechanisms 
, proposed in the literature. A first step in this direction can be found in the work by Brandenburg et al, who 

^ i' studied an extension of Sandar's model including bias but no spatial dependence (2^ . By intrinsic we mean either 
1^ [ biasing effects related to fundamental physics, such as parity- violating weak neutral currents (WNC) (TgI . [35j or - 
; I ■ given that we know little of early-Earth's prebiotic chemistry and even less of other possible life-bearing planetary 
' platforms [25j - to some as yet unknown chemical process. By extrinsic we mean possible environmental infiuences, 
: such as circularly-polarized UV light (CPL) from, for example, active star-formation regions ^23*1 or direct seeding 
^ of chiral compounds by meteoritic bombardment j3, . The spatiotemporal dynamics of the reaction network will 
• i-H , be shown to be equivalent to a two-phase system undergoing a symmetry-breaking phase transition characterized 
' by the formation of competing domains of opposite chirality. The evolution of the domain network is sensitive to 
^ ^ Earth's early environmental history and to the magnitude of the chiral bias. Using the time-scale associated with the 
■ - - I emergence of life on Earth it is possible to obtain a lower bound on the bias. In particular, it will be shown that the 
very small bias from WNC is inefficient to generate homochiral conditions. For CPL the situation is less clear due to 
uncertainty in the nature and duration of sources, but still highly unlikely. The formalism is set up to be applicable 
to any planetary platform. 

The paper is organized as follows: In section 2 we introduce the biased polymerization model and its pertinent 
rate equations. In section 3 we describe the reduced (n — 2) model and how the net chirality can be interpreted as a 
continuous order parameter satisfying an effective potential. This allows us to introduce explicitly spatial dependence 
in the study of biased polymerization. In section 4 we describe the dynamics of homochirality using techniques from the 
theory of phase transitions with mean-field Ginzburg-Landau models in the limit of no bias. In section 5 we generalize 
our results to include a small bias, describing in detail the wall dynamics in this case and the time-scales associated 
with the development of homochirality in early-Earth. We also investigate if the onset of prebiotic homochirality on 
early-Earth could have been the result of a nucleation event. We conclude in section 6 with a summary of our results 
and an outlook to future work. 
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II. MODELING BIASED POLYMERIZATION 



Sandar's model describes how long chains of homochiral polymers may evolve from a gradual build up of chiral 
monomers [13] ■ In order to reach homochirality two processes are needed: reactions must be autocatalytic so that 
longer, chirally-pure chains may be synthesized. In addition., a mechanism for chiral amplification is also needed. This 
amplification may be achieved in a number of ways [l|, [l6l . [27j . Sandars included both enantiometric cross- inhibition 
and an enzymatic enhancement catalysed by the longest chain in the reactor pool. 

Consider a polymer with n left-handed monomers, L„. It may grow by aggregating a left-handed monomer Li or 
it may instead be inhibited by the addition of a right-handed monomer Ri to either of its ends. Writing the reaction 
rates as ks and fc/, the reaction network can be written as [30|: 



Ln + Li 


2fcs 


Ln+1 


Ln + Rl 


2ki^ 


LnRl 


Ll + LnRl 


ks^ 


Lji+iRi 


Rl + LnRl 


ki^ 


RlLnRl 



(1) 

supplemented by the four opposite reactions for right-handed polymers by interchanging L ^ R. The network includes 
a substrate S from where both left and right-handed monomers are generated. The rate at which monomers are 
generated may depend on several factors. It may be due to already existing polymers with an enzymatic enhancement 
denoted here by Cli^r-) for left(right)-handed monomers. Sandars wrote Cl = Ln and Cr = Rn, where N is the 
largest polymer in the substrate. If = 2, the case we investigate here, we can model the catalytic role of dimers 
Wattis and Coveney [HI proposed instead Cl =J2-^n and Cr = J2Rn, while Brandenburg et al. [24] suggested 
a weighted sum, Cl — "-^n and Cr = ^ ni?„. Motivated by mathematical simplicity and by Soai's reactions, we 
will follow Sandars. 

Another factor that may influence the production rate of monomers is an explicit bias towards a specific handedness. 
We assign a chiral-specific reaction rate kL{R) such that the generation of left and right-handed monomers from the 
substrate S is written as 

S ^ Ll 

S ^ Rl. (2) 

The reaction rates are related to the equilibrium population of each handedness as fc^ cx exp[— E'^/fesT] and kn cx 
exTp[— {El + Ef)/kBT], where fcs is Boltzmann's constant, T is the temperature, and Ef denotes the energy difference 
between the two enantiometers, here chosen arbitrarily to suppress the right-handed monomers. Note that similar 
results would have been obtained by considering El(b.) to be the activation energy for forming L(R) molecules, and 
Ef the difference in activation energy between enantiometers. 

Kondepudi and Nelson used a similar parameterization to express the bias due to parity violation in the weak 
nuclear interactions, which has been estimated to be g = Ef/ksT ^ ^0"^^"^^ at room temperature [1, [13, H^. On 
the other hand, CPL biasing depends on a number of unknowns such as the nature of the UV source, its distance 
and duration and it's harder to estimate [1, [l3- Since in general g <C 1, one obtains kL/kji ~ 1 + g. Introducing the 
average reaction rate fee = ^^^i^, we can express the left and right-handed reaction rates as 

kL = kc{l+g/2) + 0{g^) 

kR = kc{l-9/2) + Oig^) . (3) 
If the concentration of the substrate [S] is maintained by a source Q, it will obey the equation, 

^ - - (Ql + Qr) , (4) 
where Ql(r) are the sources for left(right)-handed monomers. From eqs. [5] and [31 the sources can be written as 

Ql = kc{l + 9l2)[S]{pCL + qCR) 

Qr ^ kc{l-g/2)[S]{pCR + qCL) , (5) 



where we introduced the fidelity / of the enzymatic reactions, written in terms of p and q as p = {1 + f)/2 and 
g = (1 — /)/2. In the absence of bias {g = 0), the fidelity / controls the evolution of the net enantiometric excess. 
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As has been thoroughly discussed in the hterature [la, [2J, [30, [Ij] , the dynamical system defined by the polymer- 
ization equations for a given value of N shows a bifurcation behavior at a certain critical value of /, fc- The specific 
value of fc depends on the choice made for the enzymatic enhancements CL(^fi) and on the ratio of reaction rates 
ki/ks, but the behavior is qualitatively the same. Brandenburg et al., with kj/ks — 1, obtained fc — 0.38, with fc 
increasing with weaker cross-inhibition. In the limit /c/ — *■ 0, /c — ^ 1 and no enantiometric excess develops |24| . 

The full set of reaction rate equations governing the behavior of an n-polymer (n — 1, . . . , N) system consists of 
the following equations [23,[33|: the equation for the substrate concentration, eq. 21 the equations for the left and 
right monomers. 



d[Li] 

dt 
d[Ri] 

dt 



= Ol-Al[Li] , 

= Qr-Xr[Ri] , (6) 



where 

Al = 2ks J2 [-^«] + 2fc/ [^"1 



JV-l Af-l 



n—1 n—1 
N-1 N-1 



ks J2 i^nRl] +klY, [RnLl] , (7) 



n=2 n=2 



and 

Xr = 2ks I-R"] + t-^" 



Af-l Af-l 



n—1 n—1 
Af-l Af-l 



ks Y [-^"^i] + ''iY ' 



n=2 n=2 



and, finally, the rate eqns. for n > 2, 

d[L„ 



dt 



2ks[Li] ([i„_i] - [L„]) - 2ki[Ln][Ri] , (9) 



supplemented by the ones obtained substituting L ^ R. The factors of 2 on the rhs reflect that monomers may attach 
to either end of the chain. For n = 2, however, one must discount this for the interaction of 2 single monomers. 



III. REDUCED BIASED MODEL 

Previous authors [U, [13, [l^l have explored the evolution of the reaction network for different values of N . Here, 
we are interested in investigating not only the temporal evolution of the various concentrations ([in], [Rn]) but also 
their spatial behavior in the presence of bias. We are motivated by interesting work by Saito and Huyga '2^, who 
investigated spatial proliferation of left and right-handed polymers in the context of equilibrium Monte-Carlo methods 
and, in particular, by that of Brandenburg and Multamaki [g| (henceforth BM), where the spatiotemporal evolution 
of left and right-handed reaction networks was investigated in the absence of chiral bias. 



A. Biased Polymerization Equations 



As remarked by Gayathri and Rao [lO|, taking the concentrations to be functions of position implies that the 
number of molecules per unit volume is assumed to be large enough so that the concentrations vary smoothly with 
space and time. In other words, the concentrations are defined in a coarse-grained volume which, of course, must 
be larger than the smallest relevant distance-scale (^), to be derived below. The chemical mixture is then defined in 
block- volumes which are multiples of ~ This is essentially the procedure adopted in studying mean-field models 
of phase transitions in the Ising universality class as, for example, in ferromagnetic phase transitions, where the order 
parameter is the coarse-grained magnetization over a block of spins [2l| . Indeed, in a recent work [T^ Gleiser and 
Thorarinson (GT) demonstrated that chiral symmetry breaking in the context of the continuous model of BM can be 
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understood in terms of a second-order phase transition with a critical "temperature" determined by the strength of 
the couphng between the reaction network and the external environment. One cannot speak consistently of symmetry 
breaking without including spatial dependence. 

Adding spatial dependence to the reaction network greatly complicates its study, as we must investigate coupled 
PDEs as opposed to ODEs. Fortunately, as remarked by BM, it is possible to truncate the system to iV = 2 and still 
capture its essential behavior, the dynamics leading (or not) to homochirality within a large volume V 3> Given 
the catalytic role of dimers in the Soai reaction we will investigate the biased spatiotemporal dynamics of reaction 
networks with N = 2. 

The great practical advantage of the truncation is that it elegantly reduces the system to an effective scalar field 
theory, where the field - the order parameter - determines the net chirality in a given volume Q . 

The reaction network is further simplified by assuming that the rate of change of [L2] smd [i?2] is much slower than 
that of [Li] and [Ri\- The same for the substrate [5'], so that d[S]/dt = Q — {Ql + Qb.) — 0- This approximation 
is known as the adiabatic elimination of rapidly adjusting variables [15] : the typical time-scale for changes in the 
concentrations of larger chains such as [L2] and [R2], which depend on (or are enslaved by) the concentrations of [Li] 
and are slower than those for the concentrations [Li] and [Ri] themselves. This approximation breaks down in 
the unlikely situation where the syntheses of dimers and higher chains have similar time-scales as those for monomers. 

Using the adiabatic approximation in eq. [9] for n = 2, we can express the concentrations [L2] and [R2] in terms of 
the concentrations for the monomers as 



Also, using eq. IHwith Cl 



[L2] and Cr 



r;[^2] = 



[Ri? 



dt 

d[Ri] 
dt 



Q- 



Q- 



-2ks/ki[Riy' ^' 2[Ri] + 2ks/ki[Li] 
[R2], the equations for the monomers can be written as 



(10) 
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(11) 



Since g <C 1, we can expand the rhs of eqs. [TT] and greatly simplify the equations for the monomers. Once we 
substitute eqns. [TO] into eqns. [11] we are left with two equations for the two unknowns, [Li] and [Ri]- 

It is useful at this point to introduce the dimensionless symmetric and asymmetric variables, S = X + Y and 
A = X ~ Y , where X = [Li]{2ks/Qy^^ , and Y = [Ri]{2ks/Qy^^ are dimensionless concentrations. We can thus 
add and subtract the equations for the monomers in order to obtain the equations satisfied by the variables S and A, 
respectively. After some algebra we get, for ks/ki — 1, the biased polymerization equations. 



_idS 
" 'dt 
idA 
'dt 
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2/ 



SA 



A^ 



SA 



SA 



A^ 



(12) 



where the parameter Aq = {2ksQY^'^ , has dimension of inverse time. 5 = 1 is a fixed point: the system will tend 
towards this value at time-scales of order A^^, independently of g. With 5 = 1 and (7 = 0, the equation for the chiral 
asymmetry has fixed points at ^ = 0, ± ^/2f — 1. An enantiometric excess is only possible for / > /c = 1/2. 

Setting 5 = 1, the equation for A can be written as Aq"^^ = —dV/dA, where the dot denotes time derivative. The 
"potential" V controlling the evolution of A is an asymmetric double-well, (symmetric for g = 0) 



V{A) 



2 



f\n{l+A^ 



[A-gf 



A 



l + A" 



arctan(yl) 



(13) 



In figure 1 we show the potential V^(yl) for various values of the fidelity / and the asymmetry g. Note that complete 
chiral separation occurs only for / = 1 and is forbidden for / > 1/2 (left). The presence of asymmetry (g 7^ 0, right) 
clearly biases one of the chiralities, leaving the minima at ^ = ±1 unchanged for / = 1. 



B. Introducing Spatial Dependence 



In order to introduce spatial dependence for the concentrations, we follow the usual procedure in the phenomeno- 
logical treatment of phase transitions [l3,[2l|, by rewriting the total time derivatives in eqs. [TO]as d/dt — > d/dt—k'^^, 



FIG. 1: Left: Potential V{A) for g = and varying fidelity /. From bottom to top, / = 1,0.75, 0.5. Right: Potential for / = 1 
and varying asymmetry g. The bold line corresponds to the symmetric case, g = 0. The dash line to g = 0.1 and the dot-dash 
line to g = 0.2. 

where k is the diffusion constant. Some illustrative values of k are: k = lO^^m^s"-'^ for molecular diffusion in water 
and k — 10~^m^s~^ for air. The diffusion time-scale in a length L is r^iff — L'^/k. 

It is convenient to introduce the dimensionless time and space variables, to = Xot and xq = x-^/Ao/fc, respectively. 
The equations are then solved in terms of the dimensionless variables. Dimensionful values are obtained for particular 
choices of the parameters ks, Q, and k. For example, using as nominal values ks ^ lO^^^cm'^s"^, Q ~ 10^^cm~^s~^ 
and k for water, one obtains, t ~ 2.3 x 10~^tQ y and x ~ 8.5 x 10~'^a;om, while Aq — \/2 x 10^^s~^. 

The main consequence of introducing spatial dependence is that now the net chiral asymmetry will evolve in both 
space and time. With a racemic or near-racemic initial distribution, a typical spatial volume V ^ will coarsen 
into domains of left and right-handed polymers, separated by an interface or domain wall with approximate thickness 
^yk/Xa: an initially racemic solution gradually separates into chiral domains. The "wall" between homochiral domains 
is to be interpreted as the region of space which remains racemic [yl(x, t) = 0.] In order for an initially racemic mix 
to evolve toward homochirality, the walls separating the left and right domains must move. We will discuss in detail 
below what physical processes may trigger the wall mobility. 

In one spatial dimension, and for / = 1 and g = 0, the solution is well approximated by a "kink" profile obtained 
in the static limit {A = 0) as Akix) = tanh(— axi/ Ao/fc), where a can be found numerically to be a ~ 0.58 0. 

Using the dimensionless variables defined above, the energy of a static spatially-extended configuration in d spatial 
dimensions is given by E[A] = (fc/Ao)'*''^ / d'^Xo [^VqA ■ VqA + V{A)] , where appropriate boundary conditions must 



Gleiser and Thorarinson coupled the net chirality .4(x, t) to an external environment modeled by a stochastic 
force (C(i,x)) with zero mean {{(,) — 0) and two-point correlation function (C(x', t')C(x, <)) = a^(5(x' — x.)6{t' — t), 
where measures the strength of the environmental influence. For example, in Brownian motion, = 2jkBT, 
where 7 is a viscosity coefficient. With the dimensionless variables introduced above, the noise amplitude scales 
as Oq A(7^(Ao/fc)^/^a^. Notice that by writing 7 = A070 and identifying the noise amplitudes, 2(Ao7o)(fcBT) = 
Ao(fc/Ao)''^^a§, we obtain that the thermal energy, fc^T, has dimensions of (fc/Ao)'*^^ and thus of [length]'' as it should. 
The equations of motion in the presence of noise and no bias are 



be imposed. 



IV. CHIRAL SELECTION AS A PHASE TRANSITION I: NO BIAS 




(14) 



(15) 
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FIG. 2: Phase diagram for the 3d volume-averaged net chirality (.4)3d as a function of external noise Oenv The error bars 
denote ensemble averaging over 20 runs. 

A. Critical Point for Chiral Symmetry Breaking 

As shown by GT, in the absence of chiral bias and with / = 1 and 5 = 1 (an attractor in phase space), the 
system described by eq. [14] has critical behavior controlled by the noise amplitude a^. For a above a critical value, 
Oc, (A) — > and the chiral symmetry is restored. The brackets denote spatial averaging, (A) = (l/V) J A(Px, where 
V is the volume. In analogy with ferromagnets, where above a critical temperature the net magnetization is zero, 
one may say that above ac the stochastic forcing due to the external environment overwhelms any local excess of 
L over R enantiometers: racemization is achieved at large scales and chiral symmetry is restored throughout space. 
Thus, the history of chirality on Earth and on any other planetary platform is inextricably enmeshed with its early 
environmental history. Although in the present work we will only consider the case of "gentle" (diffusive) evolution 
of the reaction network, it should be noted that "violent" disturbances may greatly affect the results and, thus, the 
final net chiral excess. Preliminary results for 2d turbulent stirring have been recently presented, showing that indeed 
it may accelerate the emergence of a final homochiral state [61]. We are presently pursuing an alternative approach 
wherein the interactions with the environment are modelled stochastically (cf. eqs. [14] above) and hope to report on 
our results shortly. 

The equation dictating the evolution of the enantiometric excess A was solved with a finite-difference method in a 
1024^ grid and a 100^ grid with St = 10^^ and 6x = 0.2, and periodic boundary conditions. In 2d, this corresponds 
to simulating a shallow pool with linear dimensions of £ 200cm. In order to obtain the value of Uc, the system was 
prepared initially in a homochiral phase chosen to be {A){t = 0) = 1. The equation was then solved for different 
values of the external noise amplitude, a. As shown in GT, for > ~ 0.65(fc/Ao)"^^^, (^1) — > 0, that is, the 
system becomes racemized. (A) approaches a constant for large times, indicating that the reaction network reaches 
equilibrium with the environment. For d = 2, ~ 1.15(fc/Ao). In figure 2 we show the phase diagram for chiral 
symmetry restoration in 3d. 

B. Ginzburg Criterion in the Absence of Chiral Bias 

It is possible to gain much insight into the relationship between the microscopic and mean-field approaches using 
the Ginzburg criterion p^ . With g = and / = 1, the potential V^(-4) is a symmetric double- well with minima at 
A± = ±1 and a maximum at ^ = 0. A = A± are fixed points: if the system is prepared away from them, it will 
evolve towards them. When domains of both phases are present, they will be separated by domain walls. In the 
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absence of bias {g = 0) and of environmental coupling (a = 0), the only force on the walls comes from surface tension: 
the walls will straighten in order to minimize their radii within a given volume. When the average wall radius, R{t), is 
comparable to the linear dimensions of the confining volume {£) , the wall will stall or move exceedingly slowly. Given 
that left and right-handed life forms cannot coexist in the same domain, this simple model predicts, quite reasonably, 
that other factors in early-Earth's history have intervened to promote the observed homochirality. 

Within the continuous description, the smallest volume is the correlation volume, ~ 4^^, where the correlation 
length is related to the potential V{A) by = V"{A = ±1) [l^. The energy Eg required to flip a correlation- 
volume cell of a given chirality into one of opposite chirality is given by the energy barrier (AY) times the correlation 
volume: Eq = V^AV, where AV^ = \V{0) — V{ztl)\. On the other hand, if Ef is the energy to flip the chirality of a 
single molecule and is the average number of molecules in a correlation volume, then Eq = N^Ef. Equating both 

expressions for Eq we obtain, Ej — ^AV. 

From the expression for the potential, eq. [131 we obtain, with / = 1 and <? = 0, AV = 0.193 and ^ = (fc/Ao)^/^, so 
that ~ 4(A:/Ao)'^/^. In order to estimate N^, note that the parameters Q and ks define the microscopic length-scale 
Cmicro = {Q/ksY^^^- Using Q = lO^^cm^^s"^ and kg = 10~^^cm^s~^, we obtain, ^micro — 2.154 x 10"'^ cm and 
~ (C/6nicro)^ Thus, E^ ~ ^L^.^V ~ 7.7 X 10-2icm3. 

This energy can be compared with the results from the numerical study of chiral symmetry breaking in GT, who 
found ~ 0.65(fc/Ao)'^^^ — 0.4cm'^. This is the critical energy for a correlation-volume cell to flip chirality. So, we 
must divide it by the number of molecules in a correlation volume in order to obtain the critical energy per molecule, 
^num ^ g g ^ 10~^^cm^. The ratio of the two energies is, ~ 1.18, a nice agreement between the theoretical 

prediction from the Ginzburg criterion and the numerical results. 

V. CHIRAL SELECTION AS A PHASE TRANSITION II: INCLUDING BIAS 

Environmental effects, if above a certain threshold, may destroy any net chirality, restoring the system to a racemic 
state. Once external perturbations cease, the system will relax to its thermodynamically preferred state as it evolves 
toward final equilibrium. Within the present mean-field model, this evolution is characterized by a competition 
between left and right-handed domains separated by interfaces. This evolution is determined by the initial domain 
distribution and by the forces acting on the interfaces. 

A. Percolation Constraints on Chiral Bias 

One may think of the domains of each chirality as two competing populations immersed in an environment at 
"temperature" T. We use quotes to stress that this external infiuence may be attributed to several different sources 
of white noise with Gaussian amplitude . Consider a large volume V = Vl + Vr ^ Vj, where Vl{r) is the total 
volume in left(right) domains. If the fractional volumes of both left (Vl/V) and right (Vfj/V) domains exceed a 
critical value Pc, they both percolate [13] and the total volume will be a convoluted structure similar to that of a 
sponge, with regions of left and right chirality separated by a thin interface. If, instead, only one of the two fractional 
volumes exceeds pc^ the volume will be percolated by the dominant handedness. The isolated domains of the opposite 
handedness will shrink and disappear. 

Let Pl{r) be the probability that a randomly chosen correlation- volume cell will be left(right)-handed. In thermal 
equilibrium, the relative probabilities obey, Pl/pr = exp[— Ai^/fcsT], where AF is the free-energy difference between 
the two populations. Of course, if they are equally probable (no bias), pl/pr. — 1- Within the mean-field model, the 
free-energy difference between the two enantiometric phases is AF — AV^ , A = [^(^l) — where V{Anfi)) 

is evaluated at the potential minima. For the potential of eq. [T2]with / = 1, A ~ 2g{l — n / A) . 

For temperatures above the Ginzburg temperature (Tg), thermal fluctuations may drive the domains to flip their 
chiralities. As discussed in section 4.2, the associated energy scale is Eq = ksTc ~ V^^AV ll!i|. Given that we will be 
mostly interested in very small biases {AV ^ A), we may use the expression for AV obtained by taking g = Q. Thus, 
as the temperature drops below Tq, the chiral flipping is exponentially suppressed and the two populations are fixed 
by, ^\t<Tg = exp [—3^7] . From the results above, with / = 1, we obtain pl/pr\tg = exp[— 2^(1 — 7r/4)/0.193] ~ 
exp[-2.2245]. 

The ratio pl/pr.\tg sets the initial conditions for the subsequent dynamics of the system. The key point is whether 
one or both phases percolate. This is decided by comparing the probabilities Pl{r) with the critical percolation 
probability, pc- Now, pc depends on dimensionality and somewhat less on the shape of the lattice cell. As an 
illustration, and to be consistent with the numerical simulations in GT which were performed on 3d cubic lattices, 
we take Pc = 0.31, the result for cubic lattices in 3d (s^]. Using that pl + PR — 1, we obtain that for both phases 
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FIG. 
wall. 



3: Snapshot of evolution for g = 0, showing the two phases percolating the entire lattice and separated by a thin domain 



to percolate, g < gc = 0.36. This is an upper bound on chiral bias in order for both types of domains to coexist and 
percolate throughout the volume. Clearly, if <? > 0.36 only one domain percolates and chirality is firmly determined. 
However, unless some presently unknown effect strongly biases one enantiometer, g is most probably much smaller 
and both phases will initially percolate. 

In figure 3 we show an example where both phases percolate, with g — 0. The reader can verify that each phase 
crosses the entire lattice, while the two phases are separated by a thin interface or domain wall. This 2d simulation 
used a 1024^ lattice with periodic boundary conditions. For g > 0.36, only one of the two phases would percolate 
through the lattice: the minority phase would be constrained to form finite-volume domains that would shrink and 
disappear due to surface tension. We discuss wall dynamics next. 



B. Wall Dynamics in the Presence of Bias 



Once the flipping between phases ceases below Tq and percolation occurs, the wall network begins to evolve. Barring 
external influences, two forces will act on the walls: the curvature pressure Ps will act to straighten the walls, while 
the biasing pressure Pg = A will tend to accelerate the walls toward the unfavored phase. 

The curvature pressure can be written as Ps — a/R{t), where a — cro(fc/Ao)^/^ is the surface tension, with <to = 
J dxo [i(VoAfc)^ + V{Ak)] = J^^ dA[2y(A)]i/^ and R{t) is its average radius. Using eq. [Hwith / = 1 and g = 0, 
we obtain numerically that ctq = 0.789. [Corrections from a small g are quite negligible.] 

At Tc, both the average wall curvature and separation will be of order of the correlation length ^ = (fc/Ao)^/^. 
Thus, we can write 



2g(i - V4) 



0.544g 



(16) 



Unless g is unrealistically large, the motion is initially dominated by the surface pressure and walls will tend to 
straighten before the biasing pressure becomes effective. There are thus two possibilities: given a large confining 
volume with linear length-scale £, either pg becomes active before the walls straighten to £ ^ £) or after. In 
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both cases, once Pg does become active, the walls will move toward the unfavored phase so that the volume ~ £^ will 
eventually become homochiral. 

The key question is thus whether this converting mechanism has enough time to take place in early-Earth given what 
we know of its prebiotic history and the magnitude of biasing sources proposed so far. In order to answer this question, 
we write the average wall radius as R{t) — Rf){t)[k/ XaY/"^ so that Pg/ps{t) — 0.544g^o(^), with Ro{t — to = 0) — ^■ 
In the last expression we set the time to zero at wall formation time. Thus, for the biasing pressure to dominate, the 
average wall radius must grow to satisfy Rc > [0.544(7]"^: small values of g imply in very large radii. We must next 
estimate the time (tg) it takes for this critical value to be achieved. The equation controlling the wall radius is 



^ dR g 
Ao dt \R{t) 



(17) 



The solution with -Ro(O) = 1 is 

R{t)^{k/\of^[l + 2aotof\ (18) 

Note that at formation time (to = 0), the initial radius of the walls is equivalent to the correlation length and thus to 
the initial domain size as, R{0) = £, — {k/XoY^^. One can see how the initial length-scale in the problem depends on 
the two physical parameters, the diffusion coefficient k and the inverse time-scale Aq. 
Substituting the critical value for the radius obtained above and solving for time. 



1 

2cto 



- 1 



(19) 



Using the same values as before for Q, ks and diffusivity in water, we obtain (taking g <C 1), tg ~ 4.8 x 10~^g~'^y . 
[From here on, y is short for year, My(By) for millions (billions) of years, and Bya for billions of years ago.] 

If we want life's early chirality to be decided, say, within 100 million years - sometime between the culmination of 
heavy bombardment period about 4Bya [l3| and first life about 3.5 Bya [s^, we obtain a lower bound on the biasing, 
g > 7 X 10^^(100My/tg)^/^. Values of g obtained from WNC are too small to promote homochirality on early Earth. 
The situation with CPL is not as clearcut, but it seems unlikely that such values of g can be sustained long enough 
due to the variability of possible astrophysical sources and the destruction of material by unpolarized UV light 0,13 
which might pass unfiltered through Earth's prebiotic atmosphere. A more detailed quantitative study is needed. 

For g < 10~^ it would take longer than the age of the Univere (~ 14By) before the biasing becomes active. Using tg 
in the expression for R(t) above gives an estimate for the average radius of the interface by the time biasing pressure 
takes over the dynamics: R{tg) ~ (A:/Ao)^/^(ps/pg)|TG — 2(jr~^cm. 

Once biasing pressure takes over, the interface will move toward the unfavored phase. We can estimate its velocity 
by studying the equation governing its motion [2l!| . Starting with the equation describing the spatiotemporal evolution 
of the chiral asymmetry, we look for a solution moving with constant velocity w, A{x — vt). The equation becomes, 

k d'^A dV V OA 

(20) 



Aq Qj. ' 

where we assumed that the interface can be approximated by the propagation of a one dimensional front. This is 
justified considering that the thickness of the interface is of order <^ and thus much smaller than its average radius 
at t > tg. Note also that upon switching space for time, the equation describes a particle moving in the presence of 
a velocity-dependent viscous force in a potential —V{A). The solution satisfying the asymptotic boundary condition 
A{x ±oo) = ±1 is determined by the interface's velocity v. Integrating eq. [20] by parts and using that the surface 
tension o- — f dx{dA/ dxY for g ^ 1, we obtain, 

AV, fkV /O'" 



where in the last identity we used the values for the potential of eq. [131 With the same fiducial values for k and 
Ao as before, v ~ 2.5(?my~^. As an illustration, for the wall to convert a distance of 1km in lOOMy, g > 4 x 10^®. 
For the walls to sweep a distance equivalent to Earth's radius, g > 2.6 x 10^^. A small bias doubly compromises the 
conversion of racemic prebiotic chemistry to homochirality: i) the time for the bias to take over, tg ~ 4.8 x 10~^g~^y 
can easily exceed the age of the Universe even for values of the bias much larger than the ones proposed thus far; and 
ii) once the bias takes over, the distance swept by the wall, dwaii(0 — 2.5.g(t/y)m, can be exceedingly small even for 
large times of order lOOMy. 
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C. Nucleation-Induced Homochirality in Prebiotic Environments 

Given that for small bias the wall motion will not set in for a very long time, it is legitimate to ask whether the 
conversion to the final homochiral phase may happen via homogeneous nucleation [flUOl- In this case, a nucleus of 
the favored phase will thermally nucleate within the unfavored phase with a rate per unit volume T(T^g) controlled 
by the Arrhenius factor 

r ~ Ao(Ao/A:)3/2 cxp[-^g(A)/A:i3T] , (22) 

where Eg{Ab) is the energy of the so-called "bounce" or critical nucleus. In the thin- wall approximation, valid in the 
limit of small asymmetry between the two cases and thus here {A/AV <^ 1), the energy of the bounce of radius R is 
well-approximated by E{R) = —^R^AV + inR'^a. Extremizing this expression, we find, for the radius of the critical 
nucleus and its associated energy, R^ = 2a/ AV ~ lMg~^{k/Xo)^/'^ and E{Rc) = ^jfvj^ - 44.55~2(fc/Ao)3/2^ 

where in the last expressions we used the potential of eq. 1131 Note how E(Rc) oc g~^. From the expression for the 
nucleation rate we can estimate the time-scale for nucleation in a given prebiotic volume Vpb as 

Tnuc ^ (r^pb)-^ ^ 1.3 X 10-15 [-^F^J exp[89/32a2]y , (23) 

where the last expression was obtained using the usual fiducial values for k and Aq and the relation = 2kBT for 
the environmental fluctuation-inducing noise amplitude. 

Given that the nucleation time-scale decreases with volume, as an illustration let us consider a large volume of 
the unbiased phase, say a "shallow" cylindrical pool with volume V^b — 7r(10^m)^(100m) = tt x lO^m^. Within 
this volume, Tnuc — 4 X IQ-i^ exp[89/f;^a^]y. If we impose, realistically, that Tnuc < lOOMy, we obtain a bound on 
the critical nucleation barrier E(Rc)/kBT < 58.5. For small asymmetries, this bound cannot be satisfied. Indeed, 
in terms of the asymmetry g, and using that a typical noise amplitude capable of inducing sizeable fluctuations is 
~ 0.5 (cf. GT), we obtain g > 1.74, an unrealistically large value. For a volume with extension comparable to 
the Earth's radius, Vph — 7r(6.5 x lO^m)^(lOOm) = 1.3 x lO^^m^^, the bound becomes g > 1.53. We conclude that 
homogeneous nucleation cannot resolve the chirality issue. It remains to be seen if environmental disturbances may 
promote a faster nucleation rate [IH]. 



VI. SUMMARY AND OUTLOOK 



The spatiotemporal evolution of prebiotic homochirality was investigated in the context of an autocatalytic reaction 
network featuring enantiometric cross-inhibition catalysed by dimers and chiral bias. Domains of opposite chirality, 
separated by thin interfaces, compete for dominance. It was shown that the dynamics of the domain network is 
determined by the percolation properties of its initial distribution and subsequently by the two main forces acting 
on the interfaces, surface tension and chiral bias. Small biases of <? < 10~^ were shown to be inefficient to drive 
reasonably-sized reactor pools toward homochirality within presently-accepted time-scales for the origin of life on 
Earth. As a consequence, the present calculations indicate that WNCs cannot explain the observed homochirality of 
life's biomolecules. GPL remains a remote possibility, albeit current sources do not look promising either in magnitude 
and duration. Also, it should be noted that unpolarized UV may destroy any early enantiometric excess. 

The results obtained assume that the polymerization dynamics can be captured by truncating the reaction network 
to n < 2 and that the dynamics of dimers is enslaved by that of monomers (the adiabatic approximation). These 
approximations imply that complete chiral separation can only occur with perfect fidelity / = 1. Going beyond 
involves solving the complete network of spatiotemporal rate equations for larger values of n and for varying fidelity, 
a computer-intensive, but not impossible, task. However, given that the formation rate of higher n polymers will 
necessarily be slower, we believe that the results obtained here capture at least qualitatively the essentials of the more 
general case. It remains a challenging open question whether a simple transformation could be found to reduce the 
biased higher-n system to an effective field theory as done here for n = 2. 

What other possible sources of bias could have driven Earth's prebiotic chemistry toward homochirality? We cannot 
rule out the possibility that some unknown chemical bias satisfying the above bound might have been active. Another, 
highly unnatractive, possibility is that an unlikely large statistical fluctuation towards one enantiometer did occur and 
established the correct initial conditions. Possible bombardment from meteors contaminated with chiral compounds 
could also have jump-started the process, although one still needs to explain how the chiral excess formed in the 
meteors in the first place. 

It seems to us that the answer to this enigma will be found in the coupling of the reaction network to the environ- 
ment. We note again that the results obtained here are within the diffusive, and hence "gentle," evolution towards 
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homochirality. Early-Earth, however, was a dramatic environment. Given the nonlinear properties of the spatiotempo- 
ral equations describing the evolution towards homochirality, environmental disturbances, such as meteoritic impacts 
or volcanic eruptions, must have played a key role in early-Earth's prebiotic chemistry. These disturbances, if violent 
enough, would certainly affect the evolution of the chiral domain network and possibly change the bounds obtained in 
the present work. Gleiser and Thorarinson proposed to model the coupling to an external disturbance stochastically 
[l^ . Within their framework, results will depend on how the amplitude of the external "noise" compares with the 
critical value described in section 4.1. Preliminary results indicate that large enough noises (modelling external influ- 
ences) may redirect the direction of homochirality entirely, erasing any previous evolution toward either handedness. 
Further work along these lines is underway. In a different approach, Brandenburg and Multamaki suggested that 
hydrodynamic turbulence could have sped up the march toward homochirality [oj. In either case, it is clear that the 
evolution toward homochirality, as that of life itself, cannot be separated from Earth's early environmental history. 

The author thanks Gustav Arrhenius, Jeffrey Bada, Freeman Dyson, Leslie Orgel, and Joel Thorarinson for stimu- 
lating discussions. He also thanks Joel Thorarinson for producing figure 3. 
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